function tabel = multipleruns(data,num_s,num_c,num_runnings)

    together=zeros(num_s,num_s);
    for i=1:num_runnings
        cl=cluster(data,num_s,num_c);
        for j=1:num_c
            vec=cl(:,2)==j;
            vec=vec.*cl(:,1);
            vec([vec==0])=[];
            if sum(vec) > 0
                for A=1:length(vec)
                    for B=1:length(vec)
                        together(vec(A),vec(B))=together(vec(A),vec(B))+1;
                    end
                end
                vec';
            end
        end
    end

    together = together(any(together,2),any(together,1));
    %together = round(together/num_runnings*10);

    tabel=[cl(:,1) together]; tabel=[num_runnings cl(:,1)' ; tabel];

end


